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Electronic properties of heterostructures in which a finite number of Mott-insulator layers are 
sandwiched by semi-infinite metallic leads are investigated by using the dynamical-mean-field 
method combined with the Keldysh Green's function technique to account for the finite bias voltage 
between the leads. Current across the junction is computed as a function of bias voltage. Electron 
spectral functions in the interacting region are shown to evolve by an applied bias voltage. This 
effect is measurable by photoemission spectroscopy and scanning tunneling microscopy. Further 
predictions are made for the optical conductivity under a bias voltage as a possible tool to detect a 
deformed density of states. A general discussion of correlated-electron based heterostructures and 
future prospect is given. 
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I. INTRODUCTION 

Correlated-electron materials including transition- 
metal oxides have provided the basis for a variety of 
properties such as high-T c superconductivity in cuprates 
and colossal-magnetoresistance (CMR) in manganites. 1 
Therefore it is anticipated that these materials will be- 
come the core of future electronic devices. For this pur- 
pose, the fabrication and characterization of interfaces 
including correlated-electron materials are of fundamen- 
tal importance. In recent years, we have seen a tremen- 
dous amount of work on such heterostructures. Develop- 
ments in thin film growth technique- enable us to fab- 
ricate "digital" heterostructures out of correlated mate- 
rials with atomic resolutions^ and developments in 
scanning microscopy- 7 ^! 9 - enable us to characterize the 
physical properties including conduction band electron 
distribution of the heterostructures^ Furthermore, fab- 
rication of practical devices has already started. This in- 
cludes Josephson junctio n 10 : 11 and tunneling magnetore- 
sistance (TMR) junction s 12 : 13 with device characteristics 
that remain to be optimized. 

These experimental developments stimulated a great 
deal of theoretical research on correlated-electron 
heterostructures! 14 i 15 ' 16 i 17 i 18 i 19 So far, most work has fo- 
cused on the ground-state properties such as spectral 
function, charge density distribution, and linear con- 
ductance. As metal-semiconductor (or band insulator)- 
metal heterostructures provide one of the fundamental 
building brocks of current electronics ) 20 ' 21 : 22 establish- 
ing the electric properties of correlated heterostructures 
is of great importance. Detailed analysis of the trans- 
port properties of correlated heterostructures including 
current-voltage characteristics is directly relevant to de- 
vice applications ! 10 : 11 ' 12 : 13 However, this area remains 
largely unexplored. 

In this paper, we investigate the properties of hetero- 
j unctions consisting of a correlated insulator embedded 
in semi-infinite metallic leads. In particular, we focus 
on the steady state with an applied bias voltage between 



the leads. The effect of strong correlations is treated by 
using the dynamical-mean-field theory (DMFT)22 com- 
bined with the Keldysh Green's function technique 2 ^ to 
account for the finite bias voltage. We present cur- 
rent vs voltage characteristics of the junctions, position- 
dependent spectral functions, and optical conductivity 
spectra, and discuss how the transport properties of het- 
erostructures are affected by correlations. 

The paper is organized as follows: In Sec. |TT1 the the- 
oretical model, formalism, and numerical techniques are 
outlined. In Sec. IIII1 we present numerical results, and 
finally in Sec. IIV1 we discuss related work and future 
prospects. 



II. MODEL AND FORMALISM 



Model 



We consider electrons moving on a cubic lattice with 
lattice constant a (=1) and discrete translational in- 
variance in the xy plane. Thus each site is labeled by 
r = {r\\,z) with rj| = (x,y). A Hubbard-type interac- 
tion U is introduced at a number N of layers (sample S) 
located from z = 1 to N, and noninteracting leads are 
located at z > N + 1 (lead R) and z < (lead L). We 
consider the nearest- neighbor transfer t (t a ) of electrons 
in the sample (lead a) , the hybridization v a between the 
sample and lead a, and the layer-dependent potential Eq 
(see Fig. [IJ. Thus the Hamiltonian for this system is 



written as H = Hs + J2 a = 
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the in-plane momentum ku and z-axis coordinate repre- 
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FIG. 1: Schematic view of the model heterostructure. (a) 
Projection on the [100] plane. System consists of a cubic 
lattice with discrete translational invariance in the xy plane. 
Sample with a finite interaction U (1 < z < N) couples to 
noninteracting leads R (z > N + 1) and L {z < 0). Trans- 
fer intensity in the sample and lead R(L) are denoted by t 
and tjitLf, respectively. Hybridization strength between the 
sample and lead R(L) is vml). (b) Potential profile across 
the junction. Note that potentials shown are Eo(z) shifted by 
U(z)/2. 



and 



H a = -t a ^2 ( d la d r'cr + H.c.) + ^2e {r)np cr , (2) 

(r,r') r,<r 

H s -R = -v R J2(di +ma dp {l+{N+1)M<T + R.c), (3) 



H S -l =~vlY^ (4,1 <r d ni +« <r + H.c.) , 



(4) 



with Hs and H a describing the interacting region and 
lead a, respectively, and H$- a the hybridization between 
the sample and the lead a. d? a is an electron annihilation 
operator at position r with spin a, and n? a = dt^d^. 
The position r in each term is constrained as explained 
above and z is the unit vector z — (0, 0, 1). 



B. Layer dynamical- mean-field theory with finite 
bias voltage 

To apply a finite bias voltage, we adiabatically turn 
on the hybridization and interaction between the leads 
and the samplej 25 i 26 Before turning on the hybridiza- 
tion and interactions, the two leads have the chemi- 
cal potentials [Mr and /it, as well as the site potentials 
£o(z > N + 1) = sr and £o(z < 0) = El- The potential 
in the central region varies accordingly. We assume that 
the two leads are infinitely large and unaffected by the in- 
teractions, and thus remain in equilibrium. Under these 
conditions, the noninteracting leads are integrated out, 
and the electronic state of the interacting region may be 
described by the following Green's function matrices in 



(6) 

Here, R(A) and K stand for retarded (advanced) and 
Keldysh components of the Green's function matri- 
ces, respectively, and G A (k^uj) — {G R (k\\,uj)} . g R 
is the noninteracting retarded Green's function given 
by = oj + i0+ - H s (k h z,z';U = 0). 

t R (k h u) and t K {k h u) are the retarded and the 
Keldysh self-energies, respectively, representing the ef- 
fects of both electron correlation and the hybridization 
with the leads. The latter part of the self-energy can 
be obtained from the noninteracting Green's functions. 
Thus what must be determined is the self-energy due to 
correlations. 

In order to fix the correlation part of the self-energy, 
we generalize the layer DMFTi 16 i 27 i 28 Dynamical-mean- 
field theory using the Keldysh technique has been applied 
to solve the bulk correlated-electron model influenced by 
a time-dependent external field , 29 ' 30 while our focus is 
the steady state. In the layer DMFT, the self-energy due 
to correlations is approximated to be diagonal in layer 
index and independent of in-plane momentum. Thus the 
lattice self-energy is written as 

E 2, Z 'C% W ) => { S IM + v R9R,(k\\,oS)5 Z}N 

+v 2 L gl(k\\,Lo)5 Ztl }, (7) 



where 7 = R, A and K . g. 



R(A), 



■h 



j) is the retarded 



(advanced) Green's function of lead a = R,L projected 
on the layers adjacent to the interacting region, and 

gZ(k h u)={l-2f a (oj)}{g R (k h u)-g A (k {l ,u>)} (8) 

is the Keldysh Green's function of lead a. This function 
describes the distribution of electrons in terms of the 
Fermi distribution function f a (u>) — {exp/3(cj — jj, a ) + 
l}- 1 with inverse temperature (3 = 1/T and chemical 
potential fj, a . 

In DMFT, the quantum impurity model is intro- 
duced as a mathematical tool to compute the electron 
self-energy. The self-consistency condition of DMFT 
is closed by identifying the impurity Green's function 
Gimp.z(<+>) with the local part of the lattice Green's func- 
tion Gi OCiZ (lo) as 



imp,z \ 



(27T) 



(dfc||) 3 GZ,(fc|,,w), (9) 



where 7 = R(A) and K. The impurity model at each 
z is characterized by the hybridization function A z (uj) 
and the nonequilibrium distribution function of electrons 
f e ff tZ (u)) (in equilibrium, this is just the Fermi distribu- 
tion function) . The hybridization function represents the 
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intersite virtual transfer of electrons in the form of the ef- 
fective conduction band coupled to the impurity orbital. 
This includes intralayer and interlayer couplings and the 
hybridization between the sample and the leads. As in 
the equilibrium case* 2 ^ the self-consistency condition de- 
termining this function is 

Af-^H = w - s (z) E?W( W ) - {g££>( W )}T(10) 

The nonequilibrium distribution function is fixed by the 
local Keldysh Green's function as 

G£o» - {1 - 2f eff , z (u>)}{Gt P Au) - G? mp Ju)}. 

(11) 

Note that in the steady state we are focusing on, there 
is no net charge flow between the impurity orbital and 
the effective conduction band, i.e., the impurity model 
is in local "equilibrium" described by f e ff,z- There- 
fore, the Keldysh components of the self-energy and the 
hybridization function are related to the retarded ones 
by Ef H - {1 - 2/ e// , z (u;)}lm£*M and Af ( w ) = 
{l — 2/ e // iZ (o;)}lmA^(a;), respectively. Thus the self- 
consistency condition of DMFT is closed by Eqs. fTO]) 
and HTJ with Eq. ©. 

The remaining task is to solve the quantum impurity 
model defined by the hybridization function A z (lo) and 
the distribution function f e ff,z(w) with the local interac- 
tion. For this purpose, we employ the equation-of-motion 
decoupling (EOM) scheme ) 31 ' 32 which has been applied 
to study the nonequilibrium properties of quantum dots 33 
and to solve the quantum impurity models of DMFT in 
equilibrium ! 34 ! 35 ' 36 In the EOM scheme, the retarded self- 
energy is given by (including spin dependence explicitly 
as a =|, I and a = —a) 

y^R ( \ = U {" " £ "( z ) ~ A* H}(rw) - EfSjLM 
za[ ] « - eo(z) A« M - E«» - 17(1 - (rw)) ' 

(12) 

with (n za ) = -l/(k/e//,»(e)ImG^ ) , iZff (£), and 
Eg» = ~JdslmA%{e)A izcT {e) 

X (w + i0+ + s - 2e (z) - t/ + lo + i0 + - e) ' ( ' 13 ' ) 

with Ai 2CT (e) = f ef f lZrT (e) and A 2z(T (£) = l^ 6 - 

We pause to make one remark on solving the 
self-consistency equations. During the computation, 
we noticed that updating the distribution function 
directly from Eq. (fTTj) by dividing both sides by 
Gfrnp z (ui) — Gf mp z {oj) does not provide smooth con- 
vergence Instead, we found it more efficient to first 
rewrite the right hand side of Eq. pTjl as 2i{l — 

2/ e// ^(c)}|G.l^ z ( W )| 2 [lmEf H + ImAf ( w )] and then 
to update / e //, z (w) from {l - 2/ e//iZ (w)}lmAf (w). 

It may be worth pointing out the limitation of the 
EOM scheme. As discussed in detail by Gros, the EOM 



scheme does not describe the correct metallic state of the 
Hubbard model at half filling. Thus, it fails to describe 
the bulk metal- insulator transition. 34 This is because the 
EOM scheme is essentially a strong coupling expansion, 
and therefore "Kondo" physics is not fully taken into 
account. Although some amount of Kondo physics is 
included by the exchange of electrons between the im- 
purity orbital and the conduction band, it disappears in 
the particle-hole symmetric case. Therefore, it is favor- 
able to take the on-site interaction U larger than the 
critical value of the bulk Mott transition at half filling. 
When it is applied to the non-integer filling, the EOM 
scheme can reproduce the resonance-peak structure at 
the Fermi level at low temperature ! 31 ' 32 However, such 
Kondo physics is not fully taken into account. There- 
fore, results become less reliable at temperatures much 
lower than the characteristic Kondo temperature. 

Although it may not be easy, some improvement in the 
impurity solver of DMFT seems possible. For weak-to- 
intermediate interaction U, generalized iteration pertur- 
bation theory (IPT) 29 ' 37 seems realistic. But including 
higher-order perturbation processes in terms of U may 
be necessary 3 ^ because the simple second-order pertur- 
bation with respect to U does not capture the nonequi- 
librium Kondo feature when it is applied to the quantum 
dot problem^ For intermediate-to-strong coupling, the 
noncrossing approximation 4 — may be applicable. How- 
ever, this method is known to exhibit artificial structure 
in the spectral function at low temperature^ 3 - Another 
possibility is to use the recently developed continuous- 
time quantum Monte Carlo method ) 41 ' 42 which could be 
applied to a wider range of parameters. This method 
is formulated on the imaginary time axis assuming the 
system is in equilibrium. Therefore in order to apply it 
to the nonequilibrium situation, it is necessary to modify 
the method to deal with real time (or real frequency). 



C. Physical quantities 

Once self-consistency in the DMFT equations is 
achieved, one can proceed to compute physical quantities 
using the "lattice" Green's functions. Here, we first con- 
sider steady-state electric currents. In the steady state 
with a potential gradient along the z direction, the cur- 
rent is uniform in space and time. Therefore one can 
measure the current through any bond along the z di- 
rection (the steady current is conserved in the present 
formalism). Using the Keldysh Green's function matrix, 
the electric current I per spin and per unit area between 
the z and z + 1 layers is computed as 

f = I/^K + ,.fi.")-G5 +1 fi.")}. 

(14) 

Another physical quantity of interest is the optical con- 
ductivity which can provide information about the dy- 
namical properties of the system. In order to obtain 



the optical conductivity, a linear coupling between the 
current and the vector potential is first introduced ac- 
cording to the Pires-phase approximation as -I(t')A(t'). 
Then, the expectation value of the current (Ii(t)) to lin- 
ear order in Ay (t r ) is computed using the nonequilibrium 
Green's functions. Here, 1,1' — x,y,z, and c is the ve- 
locity of light. The optical conductivity tensor aiy(u>) is 
then obtained by aw(u)) — (Ii(u))) j Ey (u) where (Ii(co)) 
is a Fourier transform of (Ii(t)), and Ey(uS) is an electric 
field given by Ey(u>) — iwAy(ui)/c with Ay(u>) being a 
Fourier transform of Ay(t). When the current and/or 
vector potential is along the z direction, one has to sub- 
tract the contribution of the static current I z ^ 0. When 
both the current and the vector potential are perpendic- 
ular to the z direction, the optical conductivity is simply 
given by a Fourier transform of the current-current corre- 
lation function along the Keldysh contour. In this case, 
the conductivity a xx is expressed using the interacting 
Green's functions as^ 

z.z' 

-G-+(fc| h a;)G++(fe||,a;+ W ')},(15) 

where G = {G K + G R + G A }/2, &~+ = {G K — G R + 
G A }/2, and G++ = {G K — G R — G A }/2. In Eq. QT®, we 
neglected the vertex correction which should become sig- 
nificant in finite dimension. Nevertheless, we believe that 
the essential features of the conductivity are captured by 
the interacting Green's functions. One can use a similar 
procedure to compute other dynamical quantities such as 
the dynamical spin susceptibility. 

III. RESULTS 
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FIG. 2: Charge densities as functions of bias voltage for N = 6 
HnhharH hetprrtstmrtiirp with J J = 1 5£. 
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FIG. 3: Current-voltage characteristics of Hubbard het- 
erostructures with N = 4, 6, and 8 with U = I5t. Note that 
results shown are current per spin. Inset: current as a func- 
tion of electric field defined by E = eV/ (N + 1) . The current 
at low voltage eV/t < 14 falls onto the universal curve indi- 
cated by a light line I = ai_Eexp(— a^lE) with ai ; 2 fitting 
parameters. 



In this section, we consider the particle-hole symmetric 
case when the bias voltage is absent to ensure the Mott 
insulating state in the interacting region. The chemi- 
cal potentials and the site potentials of leads R and L 
are changed by a bias voltage V as /ii = £i = eV/2 = 
—fiR = —sr. We assume that the potential in the cen- 
tral region varies linearly as £o(z) = —Uj2 + eV(N + 1 — 
2z)/(2N + 2) interpolating el and er,. This assumption 
may be justified when the central region is highly insulat- 
ing and the density profile is not changed by an applied 
bias voltage. A more realistic calculation requires the 
potential profile to be determined self-consistently by in- 
cluding long-ranged Coulomb repulsion. This is beyond 
the scope of the current study, but some discussion will 
be given later. We mainly use the parameters U = 15t, 
t R = t L = 2.5t, v R = v L = t, and T = 0.H for the Hub- 
bard heterostructures. Qualitative behavior does not de- 
pend on the choice of parameters unless the interaction 
parameter U becomes as small as the bare bandwidth 



(and/or the interacting region becomes too thin) and the 
interacting region becomes a "metal" even at half-filling. 

Let us first check if the assumption of linear potential 
is reasonable or not. Figure plots the charge densities 
(n z ) = ^2 a {n az } as a function of applied bias voltage 
for a N = 6 heterostructure with U = 15i. As can be 
seen, the charge densities are essentially unchanged up to 
eV ~ 14i (deviation from (n z ) = 1 is smaller than 1%; it 
is about 6% even at high voltage eV^ = 20i). Therefore, 
the assumption of a z-linear potential is rather realistic. 

Next we discuss the current-voltage (I-V) characteris- 
tics of the Hubbard heterostructures. Numerical results 
are shown in Fig. [3J The thickness of the interaction re- 
gion is changed as N — 4, 6, and 8. For small voltages 
eV < 5t, current is exponentially small. The current 
grows rapidly above eV ~ 5t, and it continues up to 
eV ~ lit, above which it begins to decrease. Such a 
behavior is similar to that of a conventional metal-band- 
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FIG. 4: Layer-resolved spectral function of electrons for N = 6 Hubbard heterostructure with several choices of bias voltage 
indicated. On-site interaction is U = 15t. Shaded area shows the region occupied by electrons, i.e., f e ff tZ (u)A z (w). In (c) and 
(d), quasiparticlelike sharp structures can be seen in A z (uj) at u> = ±eV/2 corresponding to the chemical potentials of leads L 
and R. 



insulator-metal junction with the gap amplitude ~ ht 
and the distance from the middle of the gap to the top of 
the conduction band (or the bottom of the valence band) 
14i. Thus the z-axis transport is expected to be due to 
the interband Zener tunneling as discussed in Ref. [3 in 
particular in the low voltage region. In Ref. [44], break- 
down of one-dimensional Mott insulators by an applied 
voltage is investigated using a finite-size Hubbard ring. 
The applied voltage is introduced by a time-dependent 
vector potential. 

This picture is further supported by rcplotting the cur- 
rent as a function of electric field E = eV/(N + 1). The 
result is shown in the inset of Fig. [3J As can be seen, 
the current at the low voltage region (eV < 14t) falls 
onto the universal curve given by / = ai-Eexp(— a^/E) 
with 01,2 fitting parameters as indicated by a light 
line. Similar results have been reported by Al-Hassanich 
and co-workers who performed time-dependent density- 
matrix-renormalization group (DMRG) studies on one- 
dimensional heterostructures with a Mott-insulating re- 
gion sandwiched by noninteracting leads4^ This may in- 
dicate that the conventional band-insulator-like trans- 



port behavior is common for the Mott insulator in all 
dimensions. However, there has not been a direct obser- 
vation confirming the Zener tunneling in correlated in- 
sulators. This is because the numerical methods used in 
the previous studies deal with the time evolution of the 
ground state under a time-dependent external field. So 
it is difficult to see the electron spectral function in the 
steady state. On the other hand, the present formalism 
directly deals with the steady state. 

Figure 0] shows the layer-resolved spectral functions 
A z (ui) = — ilmG; oc , z (cj) for N = 6 Hubbard heterostruc- 
ture with several choices of bias voltage. Other parame- 
ters are the same as in Figs. [2] and [3] Occupied regions, 
i.e., f e ff tZ A z (uj), are shaded. Figures Q] (a)- @Jc) confirm 
the naive interpretation based on the conventional band 
insulator; about 5t of Mott gap is evident from (a); a 
chemical potential of lead L (eV — 5t) is located inside 
the upper Hubbard bands (b) , and touches the top of the 
upper Hubbard band of the layer at z = 6 at eV ~ 14t 
(c). Above eV ~ 14t (d), another tunneling process sets 
in from the top of the upper Hubbard band for layers 
z < 6 to the unoccupied state of lead R, resulting in 
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FIG. 5: Optical conductivity spectra of N — 6 Hubbard het- 
erostructure with several choices of bias voltage indicated. 
On-site interaction U = 15t. 



a negative differential conductance. At low bias voltage 
eV < 14i, the deformation of spectral functions is very 
weak because the amount of electrons injected in the up- 
per Hubbard band (and holes injected in the lower Hub- 
bard band) is too small (see also the voltage dependent 
charge densities in Fig. [2]). Therefore the magnitude of 
the gap controls the z-axis transport which is less sensi- 
tive to the detail of the model. This confirms the Zener 
tunneling mechanism. 

In addition to the "rigid" shift of the layer-dependent 
spectral functions, quasiparticle-likc sharp features are 
visible in the spectral functions at u> ~ ±eV/2 in the 
high bias voltage region eV > lAt as seen in Figs. 0] 
(c) and 0] (d). In the equilibrium case, such quasipar- 
ticle features are observed only at the Fermi level. In 
the current nonequilibrium case, it is possible to have 
two "quasiparticle" peaks at u> ~ ±eV/2 — [jll,^b. a s 
in the nonequilibrium quantum doti^ The existence of 
quasiparticles may indicate metallic transport properties, 
i.e., larger current densities. However, with the param- 
eters used here, quasiparticle features appear at rather 
high voltage and the additional tunneling process sets 
in as discussed above. Therefore the current density is 
not enhanced. It would be very interesting to investigate 
the nonequilibrium behavior of correlated heterostruc- 
tures with a smaller parameter U (close to the critical 
value for the bulk Mott transition). If the "quasiparti- 
cle" structures appear before the chemical potential 
exceeds the upper edge of the spectral function at the 
rightmost layer, it would induce a larger current, result- 
ing in highly nonlinear I-V characteristics. Such a study 
requires impurity solvers suitable for weaker interactions. 

Layer-resolved spectral functions deformed by an ap- 
plied bias voltage are essentially measurable by photoc- 
mission spectroscopy (PES) and scanning tunneling mi- 
croscopy (STM). Yet, in light of the available spatial 
resolution, the latter seems plausible. Since, a current 
is already injected by an applied voltage, the insulat- 
ing nature of the sample would not be a problem for 



STM. Although indirect, optical conductivity measure- 
ments might also be useful to investigate deformed spec- 
tral functions. Using Eq. (fl~5)l . in-plane optical conduc- 
tivity spectra a xx {u>) are computed for N = 6 Hubbard 
heterostructure with several choices of bias voltage as 
shown in Fig. [5J Other parameters are the same as in 
Figs. [2HH At low bias voltage, the spectrum is domi- 
nated by an inter-Hubbard-band transition appearing at 
uj ~ U = 15t. With increasing bias voltage, the peak 
position of the inter-Hubbard-band transition shifts up- 
wards because of the change in the potential profile in 
the interacting region. The weight of the inter-Hubbard- 
band transition is reduced and transferred to the low fre- 
quency region. The Drude-like low frequency structure 
is rather broad until sharp resonancelike structures be- 
come clear in the spectral functions. In the actual optical 
conductivity measurement, one needs a rather thick in- 
teracting region compared with the radius of the incident 
light. Otherwise, the huge Drude response from the non- 
interacting leads would hide the low frequency features 
coming from the interacting region. 



Comparison with Falicov-Kimball model 

In order to see the effect of the change in spectral 
functions, in particular the evolution of quasiparticle 
peaks, on the current-voltage characteristics, we have 
also applied the current DMFT method to the spinless 
Falicov-Kimball (FK) model. In this model, itinerant c 

fermions described by — tFK S(r,r') (^f^f + H.c.) inter- 
act with localized / fermions via the on-site Coulomb 
interaction Ufk Ylf^f^ffff^ It is known that, when 
(n{) = {flfz) 7^ 0,1, the FK model does not exhibit 
a quasiparticle feature unlike the Hubbard model4£ We 
take the local interaction as Ufk = 15t and the transfer 
intensity of c fermions as tFK = l-5t so that the positions 
of the upper- and lower-Hubbard bands and the magni- 
tude of Mott gap become similar to those of the Hubbard 
model with U — 15i studied in the previous section. 

The c-fermion self-energy of the FK model has the 
same form as the electron self-energy of the Hubbard 
model T,f a {ijS) given in Eq. (fl"3j) with S 1 ( 2 ) ZC r(w) = and 
in Z a) replaced with For the computation, we fixed 

the mean density of / fermions at each layer as (n{) = 0.5 
for all applied bias voltage. 

Numerical results for the layer-resolved spectral func- 
tions of itinerant c fermions for N — 6 FK-model het- 
erostructure are shown in Fig. [5] with several choices of 
bias voltage indicated. As seen in Fig. [6] (a), the inter- 
Hubbard-band peak-to-peak distance and the gap ampli- 
tude of the FK heterostructure are almost identical to 
those of the Hubbard heterostructure. It is also evident 
that the spectral functions of the FK heterostructure are 
less sensitive to the bias voltage and do not show the res- 
onance feature unlike the Hubbard heterostructure (see 
Fig. H] for comparison). In the FK model, electric current 
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FIG. 6: Layer-resolved spectral functions of itinerant c fermions for N = 6 Falicov-Kimball heterostructure with several choices 
of bias voltage indicated. Parameters are the on-site interaction Ufk = 15t and hopping amplitude t fk = 1.5t. Shaded area 
shows the region occupied by the c fermions. For comparison, spectral functions of N = 6 Hubbard-model heterostructure with 
U = 15t and eV — are also shown as light lines in (a). Unlike Hubbard heterostructure, a quasiparticlelike sharp structure 
does not appear in A z (ui) at any bias voltage. 
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FIG. 7: Current-voltage characteristics of TV = 6 Falicov- 
Kimball heterostructure with Ufk = 15t and tFK = l-5t. 
The electric current is carried by itinerant c fermion. For 
comparison, the current per spin as a function of bias volt- 
age of N = 6 Hubbard heterostructure with U = 15t is also 
shown. 



is carried only by the itinerant c fermions. Numerical 
results for the current-voltage characteristics of N = 6 
FK heterostructure are shown as filled circles in Fig. [7] 
For comparison, the current per spin as a function of 
bias voltage for N = 6 Hubbard heterostructure with 
U = 15t is also shown as open circles. At first sight, 
the similarity in the I-V characteristics is remarkable be- 
tween the Hubbard and the FK heterostructures. This 
clearly demonstrates that the z-axis transport properties 
of metal-Mott-insulator-metal heterostructures is domi- 
nated by the inter-Hubbard-band tunneling. Yet, about 
10% (20%) enhancement of the current at eV < 15t 
(eV ~ 17 1) can be seen in the Hubbard heterostructure 
despite a factor 2/3 smaller transfer intensity than in the 
FK heterostructure. This originates from the faster nar- 
rowing of the Mott gap by an applied bias voltage in the 
Hubbard heterostructures than in the FK heterostruc- 
tures. To see this behavior semiquantitatively, let us de- 
fine the gap amplitude by the width in which the spectral 
function becomes smaller than 0.001. Note that the spec- 
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tral function is finite everywhere and, strictly speaking, 
there is no gap at finite temperature and bias voltage. 
By this definition, the gap amplitude at eV — lit is 
estimated to be 1.60t at layer z = 3 in iV = 6 Hub- 
bard heterostructure, and 2.66i at layer z = 3 in N = 6 
FK heterostructure, i.e., the former is about 40% smaller 
than the latter. Thus the effect of dynamical fluctua- 
tions in the Hubbard heterostructure is rather strong, 
overcoming the narrower bandwidth and enhancing the 
current amplitude. 



IV. SUMMARY AND DISCUSSION 

Summarizing, we have investigated the electronic 
properties of metal-Mott-insulator-metal heterostruc- 
tures under an applied bias voltage between the metal- 
lic leads by employing the dynamical-mean-field the- 
ory combined with the Keldysh Green's function tech- 
nique. We have focused on the strong coupling region. 
In this case, the current-voltage characteristics of Hub- 
bard heterostrucures were found to be quite similar to 
that of the conventional metal-semiconductor-metal het- 
erostructures. Similar current-voltage characteristics are 
also obtained by using the Falicov-Kimball model in 
which itinerant fermions do not exhibit a quasiparticlc 
feature— These findings indicate that the electron trans- 
port in correlated-insulator heterostructures is mainly 
dominated by inter-Hubbard-band tunneling. On the 
other hand, the electron spectral functions are strongly 
deformed by an applied voltage. Deformed spectral func- 
tions under a finite bias voltage are measurable by us- 
ing (spatially resolved) photoemission spectroscopy and 
scanning tunneling microscopy. Such effects may also be 
examined by optical conductivity measurements. 



Doped case 

So far, we have considered the particle-hole symmet- 
ric situation before applying a bias voltage. Here we 
consider the hole-doped situation by increasing the po- 
tential at all layers by At with the other parameters 
unchanged. As in the undoped case, we assume that 
the potential in the central region varies linearly as 
e (z) = -U/2 + At + eV{N + 1 - 2z)/(2N + 2) inter- 
polating el and sr. In a realistic situation, this assump- 
tion is not fully justified because an external electric field 
is screened by the redistribution of carriers and the po- 
tential drop occurs only near the interface regions. In 
order to describe such a screening effect, it is necessary 
to introduce a long-range Coulomb interaction and to de- 
termine the potential profile by solving the Poisson equa- 
tion seli-consistently. These effects are neglected in this 
study. Nevertheless, the qualitative behavior including 
the linear current-voltage characteristics and the change 
in quasiparticle features presented below is expected to 
be unchanged. 



Numerical results for the layer-resolved spectral func- 
tions are shown in Fig. [H] At low bias voltage, the split- 
ting of the resonance structure can be seen [Fig. [8] (b)] . 
This is similar to the nonequilibrium Kondo effect in the 
quantum dots^ With increasing V, higher peaks are 
moved inside the Mott gap, and only the lower peaks 
remain visible. Because the electric current is carried 
by the quaiparticles in the metallic region, the change 
in the spectral function by a bias voltage may affect the 
current- voltage characteristics. 

The current-voltage characteristic of the doped Hub- 
bard heterostructure is shown in Fig. [5] Near-linear volt- 
age dependence of current can be seen at low voltage 
eV < It with the change of slope at eV ~ t as indicated 
by dashed lines. This change is due to the disappearance 
of one of the two resonance peaks. 

The frequency dependent electron occupation in the 
doped Hubbard heterostructure is changed more drasti- 
cally. As shown in Fig. [5] (e), holes are injected into the 
lower Hubbard band more strongly than in the undoped 
case (see Fig. [4] for comparison). Thus the electron den- 
sity becomes completely depressed in some frequency re- 
gions (see layers at z — 3, 4, 5). These features may affect 
excitations, i.e., the optical conductivity, more strongly 
than in the undoped case. 

Numerical results for the in-plane optical conductiv- 
ity of the doped Hubbard heterostructure are shown in 
Fig. [TO] for several choices of bias voltage indicated. Con- 
trary to the Mott-insulator heterostructure, an applied 
bias voltage destroys the low frequency Drude-like peak 
efficiently. An applied bias voltage also reduces the total 
weight correlated with the reduction of electron density 
inside the interacting region as shown in Fig. 111! The 
reduced electron density strongly depends on the layer, 
indicating the emergence of a "dipole" moment. There- 
fore, in a more realistic calculation, it will be necessary 
to include the long-ranged Coulomb interaction and de- 
termine the potential profile self-consistentlyJ£ We are 
currently working on including such effects under vari- 
ous situations. The two-peak structure in the spectral 
function seen in Fig. [5] (b) is hardly observable in the 
optical conductivity. This may be because the electron 
distribution functions f e ff,z have a similar structure to 
the spectral functions, effectively reducing the contribu- 
tion from the interpeak transition. As can be seen from 
the shaded area, f e ff,z has a two-step structure. Further, 
the finite imaginary part of the self-energy broadens the 
structure in the optical conductivity even if it exists. 

Another interesting observation is the emergence of 
the negative conductivity a xx (uj) < (see the result for 
eV = 12t in Fig. [TO]). The negative conductivity should 
be remedied by including the vertex correction. Yet, it 
implies the optical response of the correlated system be- 
comes unusual in the nonequilibrium situation. This may 
be an interesting future problem. 
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FIG. 8: Layer-resolved spectral function for iV = 6 doped Hubbard heterostructure with several choices of bias voltage indicated. 
On-site interaction is U = 15i, on-site potentials at all layer are shifted by 4t while the chemical potential at equilibrium was 
fixed to 0. Thus some amount of holes is doped into the interacting region. Shaded area shows the region occupied by electrons. 
In (b), eV — t, two quasiparticlelike sharp structures are visible in A z (uj) at to = ieV/2 corresponding to the chemical potentials 
of leads L and R. On the other hand, in (c)-(e), only one structure can be seen at ui = — eV/2. The other one, which is expected 
at uj = eV/2, is moved up inside the Mott gap and therefore becomes invisible except for a steplike feature seen at z = 1 in (c). 



Other work and future prospects 

Nonequilibrium transport properties of one- 
dimensional metal-Mott-insulator-metal heterostruc- 
tures have also been studied by Yonemitsu4£ His 
main focus is on the field-effect carrier injection by 
an applied gate voltage, and he points out that the 
ambipolar current-voltage characteristics is a general 
feature of metal-Mott-insulator-metal heterostructures. 
In contrast to the behavior of conventional metal-band- 
insulator-metal heterostructures, ambipolar behavior 
is found to be insensitive to the difference in the work 
function between the Mott insulator and the metallic 
leads. On the other hand, in light of our results, Mott 
insulators and band insulators are expected to behave 
quite similarly. Thus the difference comes from the 
absence of long-range Coulomb interaction and the 
difference in the work function between the sample and 
leads in our case. Including these effects is certainly 



necessary for applying the current DMFT method to the 
more realistic system. Yet, there may appear further 
interesting phenomena in higher dimension; in this case, 
quasiparticle structure at the Fermi level emerges more 
easily than in the particle-hole symmetric case studied 
here. 

Another work on one-dimensional heterostructures was 
performed by Oka and Nagaosa»iI They applied the 
DMRG technique to study the charge density redistribu- 
tion in the interface between noninteracting metal and 
correlated-metal in the presence of work function differ- 
ence between the two. Even though they are dealing 
with one-dimensional models in which quantum effects 
are strongest, density profiles were found to be well re- 
produced by using the classical charge with the appro- 
priate gap at an integer-filling region. This observation 
about the static properties is actually consistent with the 
previous studies on the higher-dimensional heterostruc- 
tures. It has been shown that the static charge profile 
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FIG. 9: Current-voltage characteristics of N = 6 Hubbard 
heterostructure with U = 15t and on-site potential shifted by 
At. 
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FIG. 10: Optical conductivity spectra of N = 6 Hubbard 
heterostructure with several choices of bias voltage indicated. 
On-site interaction U = 15t and on-site potential shifted by 
At. 



computed by DMFT is almost identical to the one by 
the mean- field approximation. 16,28 Based on this finding, 
they proposed an interesting "interface Mott transition" 
to explain the colossal electroresistance, i.e., the large 
switching of resistance by an applied bias voltage, ob- 
served experimentally ! 49 ! 50 Their assumption is that the 
filling-controlled Mott metal-insulator transition occurs 
by applying a bias voltage, and that the Mott transition 
is of the first order in higher-than-one dimensions. The 
Mott metal-insulator transition in the higher dimension 
is characterized by a collapse of the small energy scale 
which is roughly the quasiparticle bandwidth. 23 There- 
fore, in order to have a first order transition, the tempera- 
ture must be smaller than this energy scale. On the other 
hand, applying a bias voltage V has a similar effect to 
increasing the temperature. Therefore if the bias voltage 
required to accomplish the integer filling at some layer 
becomes larger than the small energy scale, the metal- 
insulator transition would become of the second order or 
just a crossover. A very interesting theoretical question 
is whether an interface Mott transition is really possible. 

In this paper, we have focused only on a paramag- 
netic phase and "Mott" physics. However in general, 
correlated electron systems have a strong tendency to- 
ward magnetic ordering. A single-band Hubbard model 
is known to exhibit an antiferromagnetic ordering at half- 
filling. This ordering is affected by applying a magnetic 
field, changing the chemical potentials by a gate voltage, 
or injecting a current by applying a bias voltage. In any 
of these cases, the current- voltage characteristics are ex- 
pected to be modified from the results obtained in this 
paper. Therefore including magnetic symmetry breaking 
is one of the interesting extensions of this work. It is also 
possible and desirable to consider magnetic leads instead 
of non-magnetic ones. This is of the direct relevance to 
the TMR effect. 
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FIG. 11: Charge densities as functions of bias voltage for 
N — 6 Hubbard heterostructure with U = 15t and on-site 
potential shifted by At. 



Further, application of the present DMFT method 
to other models is highly desirable. This includes the 
double-exchange model (with electron-phonon coupling) 
for CMR manganites, and the multiorbital Hubbard 
model for general transition-metal oxides. For either 
model, fluctuation, ordering or melting of internal de- 
grees of freedom would affect the transport properties of 
heterostructures. Controlling these degrees of freedom 
by external fields including a bias voltage and exploring 
the new phenomena that arise are important and urgent 
tasks. 

Another possible extension is including spatial cor- 
relation beyond single-site DMFT. As the mathemati- 
cal structure of the EOM method used in this paper is 
quite similar to that of so called correlator projection 
method^ generalizing the current method to the multi- 
site problem would not be so difficult. 
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